GEODESICS IN A QUASI-SPHERICAL 
SPACETIME: A CASE OF 
GRAVITATIONAL REPULSION 

L. Herrera* 
Escuela de Fisica, Facultad de Ciencias, 
Universidad Central de Venezuela, 
Caracas, Venezuela. 



*Postal address: Apartado 80793, Caracas 1080 A, Venezuela; e-mail address 
laherrera@t elcet . net . ve 



1 



Abstract 



Geodesies are studied in one of the Weyl metrics, referred to as 
the M-Q solution. First, arguments are provided, supporting our 
beUef that this space-time is the more suitable (among the known 
solutions of the Weyl family) for discussing the properties of strong 
quasi-spherical gravitational fields. Then, the behaviour of geodesies 
is compared with the spherically symmetric situation, bringing out the 
sensitivity of the trajectories to deviations from spherical symmetry. 
Particular attention deserves the change of sign in proper radial ac- 
celeration of test particles moving radially along symmetry axis, close 
to the r = 2M surface, and related to the quadrupole moment of the 
source. 
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1 Introduction 



As we know, all physical systems in nature (without exception) are submit- 
ted to fluctuations, and of course this well established fact also applies to 
self-gravitating systems. Accordingly we have to accept that any physical 
property of such a system, e.g. spherical symmetry (if assumed ) is to be 
submitted to such fluctuations. In this paper we shall deal with two obvious 
questions arising from these general comments, namely: 

• How to describe such fluctuations?. 

• How sensitive are the properties of compact objects to deviations from 
spherical symmetry?. 

Now, if the fleld produced by such self-gravitating systems is not particu- 
larly intense (the boundary of the source is much larger than the horizon) and 
the deviation from spherical symmetry is slight, then there is not problem 
in representing the corresponding deviations from spherical symmetry (both 
inside and outside the source) as a suitable perturbation of the spherically 
symmetric exact solution [1] (although strictly speaking the term "horizon" 
refers to the spherically symmetric case, we shall use it when considering 
r = 2M surface, in the case of small deviations from sphericity). 

However, as the object becomes more and more compact, such perturba- 
tive schemes will eventually fail close to the source. Indeed, as is well known 
[2], though usually overlooked, as the boundary surface of the source ap- 
proaches the horizon (in the sense indicated above), any flnite perturbation 
of the Schwarzschild spacetime, becomes fundamentally different from any 
Weyl metric, even if the latter is characterized by parameters whose values 
are arbitrarily close to those corresponding to spherical symmetry. In other 
words, for strong gravitational flelds, no matter how small are the multipole 
moments (higher than monopole) of the source, there exists a bifurcation 
between the perturbed Schwarzschild metric and all the other Weyl metrics. 
This in turn, is just the expression of the fact [3], that the only static and 
asymptotically-flat vacuum space-time possessing a regular horizon is the 
Schwarzschild solution. For all the others Weyl exterior solutions [4], the 
physical components of the Riemann tensor exhibit singularities at r = 2M. 

In line with these arguments, let us mention that in recent works [5], 
[6] it has been shown that such bifurcation between the exactly spherically 
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symmetric case and a solution of the Weyl family, when considering the 
source (the interior), takes place for strong gravitational fields, but before 
the horizon is reached. 

It is worth mentioning that the sensitivity of the trajectories of test par- 
ticles in the 7 spacetime [7], to small changes of 7, for orbits close to 2M, 
has been brought out [8]. Also, the influence of the quadrupolc moment on 
the motion of test particles within the context of Erez-Rosen metric [9] has 
been investigated by many authors (see [10], [11] and references therein) 

Motivated by the above mentioned facts, we want to present in this work 
another approach to the problem of describing deviations from spherical sym- 
metry. This consists in considering an exact solution of Einstein equations 
(of the Weyl family, since we restrict ourselves to static axially-symmetric 
solutions) continuously linked to the Schwarzschild metric, through one of 
its parameters, instead of considering a perturbation of the Schwarzschild 
space-time. 

However, since there arc as many different (physically distinguishable 
[12]) Weyl solutions as there are different harmonic functions (see next sec- 
tion), then the obvious question arises: which among Weyl solutions is better 
entitled to describe small deviations from spherical symmetry?. 

Although it should be obvious that such question does not have not a 
unique answer (there is an infinite number of ways of being non-spherical, so 
to speak), we shall invoke a very simple criterion, emerging from Newtonian 
gravity, in order to choose our solution. 

Indeed, in order to answer the question above, let us recall [13] that most 
known Weyl solutions, present a serious drawback when describing quasi- 
spherical space-times. It consists of the fact that its multipole structure is 
such that multipole moments, higher than quadrupole, are of the same order 
as the quadrupole moment. Instead, as it is intuitively clear, the relevance 
of such multipole moments should decrease as we move from lower to higher 
moments, the quadrupole moment being the most relevant for a small de- 
parture from sphericity. Thus for example in Newtonian gravity, multipole 
moments of an ellipsoid of rotation, with homogeneous density, and axes 
(a,a,b) read: 




, — a—b 
— a ' 
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because of the factor e", this equation clearly exhibits the progressive 
decreasing of the relevance of multipole moments as n increases. 

Thus, in order to describe small departures from sphericity, by means of 
exact solutions of the Einstein equations, we would require an exact solution 
whose relativistic multipole structure shares the property mentioned above. 
Fortunately enough, such a solution exists. 

Indeed, there is one (exact) solution of the Weyl family (hereafter referred 
to as the M-Q solution [13]), which may be interpreted as a quadrupole 
correction to the Schwarzschild space-time (see below). It is for this exterior 
metric that we are going to study the motion of test particles (more exactly, 
for a sub-class of this solution). 

As we shall see, for very strong fields, arbitrarily small values of the 
quadrupole moment of the source, introduce important differences in the ra- 
dial motion of test particles, as compared with the exactly spherically sym- 
metric case. Particularly interesting is the appearance of repulsive forces, 
related to the quadrupole moment, acting on radially moving particles, along 
the symmetry axis and close to the horizon. 

The paper is organized as follows. In the next section we give a brief 
resume of Weyl solutions. Next, we describe the M-Q solution [13] and 
comment on its properties, in Section 3. The geodesic equations for the M- 
Q solution are found and analyzed in Section 4. Finally results are discussed 
in the last Section. 

2 The Weyl metrics 

Static axysymmetric solutions to Einstein equations are given by the Weyl 
metric [4] 

ds" = e^'^de - e-2*[e2r(cip2 + dz") + p'dct>\ (1) 
where metric functions have to satisfy 

*,pp + p-i*,p + *,,, = 0, (2) 

and 

V^, = p{^l-^l) ■ F,, = 2pv^,,v^,,. (3) 

Observe that (2) is just the Laplace equation for \E' (in the Euclidean 
space), and furthermore it represents the integrability condition for (3), im- 
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plying that for any "Newtonian" potential we have a specific Weyl metric, a 
well known result. 

The general solution of the Laplace equation (2) for the function pre- 
senting an asymptotically flat behaviour, results to be 

^^-t^.PnicOSe), (4) 
n=0 ' 

where r = (p^ + z^Y^"^ , cos9 = z/r are Weyl spherical coordinates and 
P„(cos^) are Legendre Polynomyals. The coefficients arbitrary real 

constants which have been named in the literature "Weyl moments", al- 
though they cannot be identified as relativistic multipole moments in spite 
of the formal similarity between expression (4) and the Newtonian potential. 

Another interesting way of writting the solution (4) was obtained by Erez- 
Roscn [9] and Quevedo [14], integrating equations (2, 3) in prolate spheroidal 
coordinates, which are defined as follows 

r+ + r- ~ ^- 

X > 1 ,-!<!/< 1, (5) 

where a is an arbitrary constant which will be identified later with the 
Schwarzschild's mass. The prolate coordinate x represents a radial coor- 
dinate, whereas the other coordinate, y represents the cosine function of the 
polar angle. 

In these prolate spheroidal coordinates, ^ takes the form 

oo 

*-5:(-ir+^g„g„(x)p„(i/) (6) 

n=0 

being Qn{y) Legendre functions of second kind and g„ a set of arbitrary con- 
stants. The corresponding expression for the function F, has been obtained 
by Quevedo [14] 

A sub-family of Weyl solutions has been obtained by Gutsunaev and 

Manko [15], [16] starting from the Schwarzchild solution as a "seed" solution. 

Both sets of coefficients, a„ and g„, characterize any Weyl metric [14]. 
Nevertheless, as mentioned before, these constants do not give us physical 
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information about the metric since they do not represent the "real" multi- 
pole moments of the source. That is not the case for the relativistic multipole 
moments firstly defined by Geroch [17], and subsequently by Hansen [18] and 
Thorne [19], which, as it is known, characterize completely and uniquely, at 
least in the neighbourhood of infinity, every asymptotycally flat and station- 
ary vacuum solution [20] , providing at the same time a physical description 
of the corresponding solution. 

An algorithm to calculate the Geroch multipole moments was developed 
by G.Fodor, C. Hoenselaers and Z. Perjes [21] (FHP). By applying such 
a method, the resulting multipole moments of the solution are expressed 
in terms of the Weyl moments. Similar results are obtained from Thome's 
definition, using harmonic coordinates. The structure of the obtained relation 
between coefficients a„ and these relativistic moments allows to express the 
Weyl moments as a combination of the Geroch relativistic moments. 



2.1 The Monopole— Quadrupole solution, M — Q 

In this section we would like to review the properties of a solution (hereafter 
referred to as the M — Q solution [13]) which is particularly suitable for the 
study of perturbations of spherical symmetry. The main argument to support 
this statement is based on the fact that previously known Weyl metrics (e.g. 
Gutsunaev-Manko [15], Manko [16], 7-metric [7], etc.) have a multipolar 
structure (in the Geroch sense) such that all the moments higher than the 
quadrupole, are of the same order as the quadrupole moment. In fact for the 
above mentioned metrics we have (odd moments are of course vanishing) 

Mf^ = _3Mf«=4?2M5 (7) 

35 

6 6 7 33 4 i5 23r ^ 7 15^^ ' 

where 52 is the quadrupole parameter in the Erez-Rosen metric. For the 
gamma metric the results are 

Mo = 7M 
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M2 


= 7^(1 




M4 


= 7^(1 




Me 


= 7— (1 





389 . 292 



495 ' 165 



For the Curzon metric [22], the situation is still worse since it posseses a 
unique parameter which represents the mass, and all higher moments are 
proportional to increasing powers of that parameter, i.e., Mq = —Cio,M2 = 

1^3 ]\J _ 19 5 ]\/r _ 389 7 

Prom the comments above it should be clear that neither of the mentioned 
metrics satisfies the criterion described in the Introduction. 

In [13] it was shown that it is possible to find a solution of the Weyl family, 
by a convenient choice of coefficients a„, such that the resulting solution 
possesses only monopole and quadrupole moments (in the Geroch sense). 
The obtained solution (M — Q) may be written as series of Gutsunayev- 
Manko [15] and Erez-Rosen [9] solutions, as follows: 

oo 

*M-Q = + + g2^g2 + . . . = ^ , (8) 

where the zeroth order corresponds to the Schwarzschild solution. 

oo \2n+l 

*.° = -E^— r^2n(cOS^) , (9) 
n=0 -T 

with A = M/ r and it appears that each power in q adds a quadrupole correc- 
tion to the spherically symmetric solution. Now, it should be observed that 
due to the linearity of Laplace equation, these corrections give rise to a series 
of exact solutions. In other words, the power series of q may be cut at any 
order, and the partial summation, up to that order, gives an exact solution 
representing a quadrupolar correction to the Schwarzschild solution. 

The simplest way to interpret physically the exact solutions obtained 
from the quadrupolar corrections described above, consists in analyzying 
the corresponding multipolar structure. Thus , it can be shown that cutting 
solution (8) at some order a, one obtains an exact solution with the following 
properties: 
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• Both, the monopole and the quadrupole moments are non- vanishing: 
Mo = M, Ma = qM^. 

• All the moments beyond M2, until order 2{a + 1) (inclusive) vanish. 

• All moments above M2(^a+i) are of order Therefore, the solu- 
tion represents a quadrupolar correction to the Schwarzschild solution, 
which is an exact solution up to the given order. 

Since we are interested in slight deviations from spherical symmetry, we 
shall consider the M-Q solution, only up to first order in q. Thus, the 
explicit solution up to first order, describing a quadrupolar correction to the 
monopole, is ( note a misprint in the equation (13) in [23]) 

= *,o + gvl>^i = iln(|^')+|g(3|/2-l)x 



X 



^ 4 37/2 - 1 J ^ Iv^TT; ~ (x2-7/2)(3^2_l) ^ ^ 



.(1) 
M-Q 



225 



gMln 



\t2 - r 

^2 _ y2 



ygx(l - 7/2) 



225 



\ 15 f o ^2 022 8 ,x2 + 1 

1 qlx^ + 7y^ - 9x^y^ + 1 - - 



32 



3 a;2 — 7/2^ 
a; — 1' 



In 



X — 1 
x + 1 



q\x' - 1)(1 - y'){x' + 7/2 - 9xV _ 1) ln2 



1024 



15 

1 -g(a;2 + 47/2 - 9^27/2 + 4) 

64 



75 
16^ 



1-y' 

2 _ 1,2 4"^'^"' ' ^ ^ ^'t.2 _ ,,2^2 



2-2 1 ^ ^g(x2 + 7/2 



x^ — y^ 



[x^ — y^ 



'^^2/06 4io42 c22|j24 4 Gn^ V'^ 

q [2x — x +6xy — ox y +4x7/ —y —y 



192 



a;2 - 7/2)4 



The first twelve Geroch moments of this solution are (odd moments van- 
ish, because of the reflection symmetry) 

Mo = M,M2 = M=^g,M4 = 0,M6 = -^MV 



(1 



9 



Ms = 

Mi2 = 



Prom the expressions above, it is apparent that the parameter q = M^/M^ 
representing the quadrupole moment, enters into the muhipole moments 
M2„, for n > 2, only at order 2 or higher. Accordingly, solution (10, 11) 
for a small value of g, up to order g, may be interpreted as the gravitational 
field outside a quasi-spherical source. The spacetime being represented by 
a quadrupole correction to the monopole (Schwarzschild) solution. This is 
in contrast with other previously mentioned solutions of Wcyl family, where 
all moments higher than the quadrupole are of the same order in q as the 
quadrupole. And therefore for small values of q they cannot be interpreted 
as a quadrupole perturbation of spherical symmetry. 

Instead of cylindrical corrdinates (p, z) , it will be useful for the next 
section to work with Erez-Rosen coordinates (r, 0) given by: 

z = (r - M) cos 9 

p = (r2-2Mr)^/2sin^ (13) 
and related to prolate coordinates, by 

M 

y = cos 6* (14) 

The metric functions of the solution, up to the first order in g, hereafter 
referred as M-Q^-*^) are: 

= ^^(l-|) + i^^3z/^-l)(3i?^-6^ + 2)ln(l-|) 

_ ^^in fl - 1) - + 1^5(3|/2 -1){R- 1) (],5) 

8^ V R) A^iR-iy-y^ 16^^^ ' ' 



M-Q 



2 



{R-lf-l 
{R-iy-y\ 



— (7 1 

8 



^,(l-y^)-^,(l-y^) 



(i?-l)M 



l)ln 



1--) 

r) 



{{R-lf-y^f 



(16) 
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with R = r/M. 

It is wortli noting that the M-Q solution (8) up to order q, imphes the 
appearance of terms of order in T^m_q as shown in (11). However such 
terms have not been included in (16) since they do not play any role in the 
discussion below. Indeed, in the case y — 1 (along the symmetry axis) all 
these terms vanish, and in the case 7/ 7^ 1, in the limit R— > 2, the leading 
term among the factors multiplying g^, is of the same order of magnitude as 
that multiplying q, therefore since we are considering small values of q, all 
terms with may be neglected. 

In the next section we shall analyze the motion of a test particle in a 
space-time described by the M-Q*^^'' solution. 

3 The geodesies for the M-Q*^^^ metric 

The equations governing the geodesies can be derived from the Lagrangian 

2jC = QapX^xf", (17) 

where the dot denotes differentiation with respect to an affine parameter s , 
which for timelike geodesies coincides with the proper time. Then, from the 
Euler- Lagrange equations, 

d ( dC\ dC _ 

dSydi^j-d^-^' ^^^^ 

it follows: 

igtt + K^9tt,r + gttfiO) = (19) 

2rgrr + 2f{fgrr,r + grrfiO) - gtt,r " r'^grr,r " 0'^gee,r " 0^5'</.</.,r = (20) 

29gee + 29{fgee,r + gee,eO) - Pgu,d - r'^grrfi - O'^geefi - <i>'^9<i><i>,e = (21) 

<P9H> + (t>iT9H,r + 9Hfi^) = (22) 
Where, our exterior metric in Erez-Rosen coordinates read: 



1-2A 



1 _ 2M^ 



g,, = -e2r-2V(l-2A + A2sin2^) = -rV(^.^-^^) 
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g^^ = -e-2Vsin2^(l-2A) = -r2sin2^e-2''V 
gtt ^ e = I 1 1 c « oi 



and the functions Fgi and "^qi are given by 
5 



^(3|/2-l)(3i?2-6i? + 2)ln(l 



2 



i?y 4 (ii- 1)2 -1/2 16 



(23) 



l^(3y^-l)(it:-l) ,(24) 



= -^(l-y^)(i?-l)ln(l-|) 



2 



(25) 



Since we are concerned only with timehke geodesies, the range of our coor- 
dinates is: 

oo>i>0 r>2M 7r>^>0 2t:> d)>Q 



3.1 Radial geodesies 

Let us consider the motion of a test particle along a radial geodesic, for an 
arbitrary value of 0. Thus putting ^ = = in (20), we obtain 



.^2/ ^ I 9rr9tt,r 9tt,r _ q 



1-2/ I tirrtm,r\ 
2rgrr + r {grr,r H ) 



9tt 



9tt 



(26) 



or using (23) 



(27) 



with ^qi, Fyi, \E'gi and Fgi ^. easily obtained from (24, 25). 

In the spherically symmetric case {q = 0) we obtain from (27) the well 
known result 

M , , 

^ = -ZJr (28) 
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Now, before proceeding farther, and in order to express our results in 
terms of physically meaningful quantitites, let us introduce (locally defined) 
coordinates associated with a locally Minkowskian observer (or alternatively 
a tetrad field associated with such Minkowskian observer). Thus, let 



and 

from where it follows that 



dX = \/—grrdr 
dT = ^tdt 



(IX 

(IT 



'dX\2 _ 1 



(29) 
(30) 

(31) 



and 



1 - 



,dX., 
-~dT' 



,dX^9 



1 - 



'dXV2 
dT ' 



dT' 2{-g„f?,/2) 



In the spherically symmetric case (g = 0), (32) reduces to 



d'^X 
dT'^ 



or, in terms of R 



fx_ 

dT^ 



M 



, dX . o 
1 - ( f 



(32) 



(33) 



1 - 



, dX . 
'~dT' 



- 9^-1/2 



{R-2) 



(34) 



which is also a known result [24]. Since (■^) is always smaller than one, the 
attractive nature of gravity for any value of r (larger than 2M) is clearly 
exhibited in (33). 

Let us now see what is the situation for any value of q, arbitrarily small 
but different from zero, as we approach the horizon. We shall analyze two 
different situations, namely: 

• The particle moves along a radial geodesies, outside the symmetry axis 



The particle moves along the symmetry axis {y^ — V) 
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3.1.1 



Let us assume that the test particle moves along a radial geodesies outside 
the symmetry axis, then as R— > 2, the leading term in (32), for small values 
of q is 



1 _ (i2L\'2 



0{q) 



(35) 



dT^ 2i(i?- 2)1/2 (i?- 2)1/2 

Observe that this limit (except for small corrections of order q) is the same 
as the one obtained in the exactly spherically symmetric situation from (34). 
Thus for particles moving along radial geodesies (excluding the symmetry 
axis), small values of the quadrupole moment introduce small perturbations 
on the trajectories, and the attractive nature of the gravitational force is 
preserved at all times. Let us now see how this situation change for particles 
moving along the symmetry axis {y^ = !)• 

3.1.2 y2 ^ 1 

In this case, it follows that as R— > 2, 



5 R,r 



-e 4(ii-2)^ 



16 (i?- 2) 

which is of course a positive quantity if g < 0, and: 

V2 



and 



5qR, 



5q 

g8(fl-2) 



5q 



2(i?-2)3 

Using (36-38) in (32), we finally find that as R- > 2 



d'^X 



-5qR, 



dT^ 2i(i?- 2)3/2 



dT' 



-5q 



g8(i{-2) 



(36) 

(37) 
(38) 

(39) 



which is a positive quantity if g < 0, diverging at the horizon, for any, differ- 
ent from zero, value of q. Thus, small deviations form spherical symmetry 
(with q < 0) produce, close to the horizon, a positive (in the sense of increas- 
ing X) radial acceleration of particles moving radially on the symmetry axis. 
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as measured by a Lorentzian observer. Therefore, under such circumstances, 
such observer would infer the existence of a "repulsive" gravitational force 
acting on the particle. This conclusion is valid in all locally Minkowskian 
frames, including of course the instantaneous rest frame where ^ = 0. Thus 
for example, if a test particle, initially at rest with respect to the Minkowskian 
local frame, is placed close to the horizon (on the axis of symmetry), as seen 
by such observer, it will tend to displace in the sense of increasing X, moving 
away from the source. 

If q > then will be negative and tending to zero as r approaches 
2M. 

It is worth noting that the source of this difference between the y"^ = 1 
and the 7^ I cases is the the term ^'gi^^- Indeed, in the later case {y"^ 7^ 1) 
we have , as R— > 2, 

15R,r(l-y^) 
^ - 16(fl-2) ' 
whereas in the case y^ = 1, we have 

5R,r 
8(R^2y 

It should be emphasized, once again, that the results above (36-39), hold 
for any small, but finite value of q. Therefore one should not expect to 
obtain the spherically symmmetric limit (34) from (39) as g — > 0, since the 
spherically symmmetric contribution (in the limit R ^ 2) being smaller, has 
been omitted in (39). 

Finally, observe that the fulfilment of equation (21), imposes a constraint 
on r which does not play any role in the very appearance of the repulsive 
effect as implied by equation (39) since, as already mentioned, such an effect is 
observed by all locally Minkowskian observers, including the instantaneously 
at rest observer. 



3.2 Circular geodesies 

Let us now consider circular geodesies. From (19)-(22), and the definition of 
the angular velocity of a test particle along circular geodesies 

. = i (42) 



*oir > (41) 
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we obtain using r = 6* = 0, in the equatorial plane {9 = 7r/2) 

In the spherically symmetric case (g = 0), (43) yields the Kepler law 



- 5 (44) 



However, for arbitrarily small (different from zero) quadrupole moment, 
(43) exhibits a bifurcation from Kepler law, as R— > 2. Nevertheless, such 
bifurcation takes place for orbits for which the tangential velocity of the 
particle is larger than 1, and therefore are deprived of physical interest (see 
below). Indeed, we have for the tangential velocity VV^ of the test particle 
along the circular geodesic (see for example [25]) 



Then, we obtain on the equatorial plane 

^2 ^ R,r+Rq'^MR'2) 

{R-2){R,r-Rq^,i,r) ^ ' 

It follows from (46) that photonic orbits (IF = 1) differ from r = 3M(i? ~ 
3) by factors of order q. Also, if R— > 2, then W'^— > Of course these 
orbits are clearly unphysical, being W > 1. 

Finally it is worth noting that circular geodesies on the equatorial plane 
are stable. Indeed, it follows from (21) with r — 9 — that 

'2Oge0^i\gu,e + u;^9^^,e) (47) 
implying 6* = on the equatorial plane. 

4 Conclusions 

We have first argued that the M-Q*^^) solution represents a good candidate 
among known Weyl solutions, to describe small deviations from spherical 
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symmetry. Then we have calculated the geodesic equations for a test particle 
moving in that space-time. 

We have found that in general the behaviour of circular geodesies is not 
very different from the spherically symmetric case. However important dif- 
ferences apppear in the behaviour of radial geodesies. Thus, it has been 
established that particles along the axis of symmetry, close to i? = 2 feel a 
repulsive force ii q < (corresponding to an oblate source). It should be 
obvious that such an effect is absent in Newtonian theory, for small {q « 1) 
values of the quadrupole moment. 

It goes without saying that we have not obtained the full description 
of the trajectories of test particles in the M-Q*^^) spacctime, since we have 
not completely integrated the geodesic equations. Certainly such integration 
would be necessary in order to extract definitive conclusions about the be- 
haviour of the test particle. Unfortunately we were not able to find analytical 
solutions to this problem, and a numerical approach is out of the scope of 
this paper. 

It is worth mentioning that repulsives forces in the context of general 
relativity have been brought out before, the best known example probably 
being the case in the Reissner-Nordstrom solution (see for example [26]). 
Also, they appear in connection with rotating and/or unbounded sources 
[27] or as spurious phenomena associated with coordinate effects (see [24] 
and references therein). 

As an example of the later is the bouncing back of particles that strike 
the surface r=2M, found by Rosen [28] for the spherically symmetric case. 
However observe that the proper radial acceleration as measured by a locally 
Minkowskian observer, for the spherically symmetric case, does not exhibit 
a change of sign in the limit r— > 2M, as indicated by (33). 

In our case, the repulsive character of the force (understood as a change 
of sign in is exhibited by means of a quantity measured by locally 

Minkowskian observers, and therefore it is not a coordinate effect, and takes 
place close to, but before reaching the surface r = 2M. In this sense the 
following comment is in order: It should be clear that at = 2 there are not 
locally Minkowskian observers, since that is a singular surface for any non- 
vanishing value of q [11]. However, for sufficiently small (but finite) values 
of q, close to but at a "finite" distance from 2M (where locally Minkowskian 
observers do exist) the leading term in is given by (39). In other words 
it is in principle possible to find a region where locally Minkowkian observers 
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exist and change of sign. 

Finally it is worth asking what are the physical implications of the study 
presented here. In view of the above, we would say that any result implying 
the presence of strong gravitational fields and spherical symmetry, should be 
tested against small fluctuations of such symmetry. Also, it remains to be 
seen which, if any, might be the implications of the reported repulsive effect, 
on some astrophysical scenarios, as for example in relativistic jets [29]. 
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